Calculation of the He-II quasiparticle spectrum by the method of collective variables 
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The method of collective variables (MCV) has been used to calculate the logarithm of the He-II 
ground-state wave function, In ^l/o, to an accuracy of a first correction to the Jastrow function and, in 
a second approximation, the wave function of the first excited state and the He-II quasiparticle 
spectrum. The functions ^'o and ^'k were found as the eigenfunctions of the A^'-particle Schrodinger 
equation, and the function ^0 was connected to the structure factor of He-II, using the Vakarchuk 
equation. The model does not contain any fitting parameter or function. The quasiparticle spectrum 
calculated numerically agrees well with the experiment. Our solution improves the result obtained 
early by Yukhnovskyi and Vakarchuk. 

PACS numbers: 67.25.dt 



I. INTRODUCTION 

The structures of the iV-particle wave functions of the 
ground and weakly excited states of helium-II are known 
in the main^'^'^'^'^'^'^'^ ' ^ ° 1 ^ ^ , and the solutions which 
take into account several first approximations have been 
obtained. In our opinion, the main unresolved problems 
concerning the microscopic physics of He-II are the struc- 
ture of the composed condensate, the nature of the A- 
transition, and the role of microscopic vortex rings. 

The form of the He-II quasiparticle spectrum has been 
forecasted by Landau for the first timc^^. In Feynman's 
known worksii^i^, an opportunity to determine this spec- 
trum making use of the structure factor has been demon- 
strated. Feynman intuitively found the structure of the 
^I^k wave function for the state of He-II with a single 
phonon and approximately deduced the He-II quasipar- 
ticle spectrum. According to Feynman and Cohen, 

^fklri, . . . ,rjv) = -iAklri, • ■ ■ ,rjv)^'o(i"i, ■ ■ • ,i"Ar), (1) 
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-kP-ki, 
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where \l/o is the wave function of the ground state, 



Pk 



1 ^ 



(k^O) 



(3) 



are collective variables^, and A'' is the total number of 
atoms in helium. However, it has not been shown in 
Feynman's works that function ^ is the eigenfunction 
of the iV-particle Schrodinger equation. Feynman's ideas 
have been developed in a great number of works (see, 
e.g.^ii^a&iiia8ii9i20)^ rpj^g Feynman-Cohen function has 

been specified in works^ where the analysis of the to- 
tal Hamiltonian of the system has been carried out. A 
more accurate form of the function V'k and the structure 
of the function have been found in works^i^iiSiiiSiii, 



where ^0 and were sought as the eigenfunctions of 
the Schrodinger equation. 

The idea of the MCV has been proposed in Bogolyubov 
and Zubarev's work"'^'^. This method has been substanti- 
ated and developed in works^i^ i^^i^"^:^^ . In workii, tak- 
ing advantage of the MCV, the ^'o and ^^k functions of 
helium-II have been calculated making use of the model 
potential of interaction between He^ atoms with one fit- 
ting parameter. Nevertheless, as was indicated in^i i^^'^^ , 
the derivation of the 5*0 and ^Pk wave functions, as well as 
the He-II quasiparticle spectrum, starting from the He-II 
structure factor, known from the experiment, rather than 
from the model potential has significant advantages. In 
this case, the problem does not contain fitting parame- 
ters, and one can avoid the task of description of atomic 
interaction at small distances, which arises because of 
atoms' extensionSiiSiii. 

Such an approach has been considered in work^, where 
In ^E'o was found in a zeroth-order approximation, while 
'0k and the He-II quasiparticle spectrum in a first one. 
The obtained spectrum E{k) agreed well with the experi- 
ment. In this work, we calculated ^"0, V'k, and E{k) more 
accurately. Namely, we found a first correction to In ^'o 
and a second ones to 0k and E{k). In doing so, we used 
the equation for '^q, derived in work^^ (below, we coin it 
as the Vakarchuk equation) . Actually, the expansion pa- 
rameter of the problem was the function 2a{k)k/ko (see 
Fig. 1), the average value of which within the interval 
fc = fco was about —1/2; i.e. the parameter was not 
small. Therefore, the corrections to In^o and E{k), gen- 
erally speaking, were not small too, and their calculation 
was of interest. 



II. THE GROUND STATE OF HELIUM-II 

A more detailed analysis of the equations and the 
method of determining ^'o were exposed in work a^i^°'^^ . 
The necessary equations for 4'o and 'ip]^ were found by 
Yukhnovskyi and Vakarchu k^i^i^°i^"^ . We shall use differ- 
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ent notations and different forms of the equations for 'i'o 
and tpyM. (the latter is partially caused by our desire to 
reduce the error of numerical solution of the equations^i) . 

The wave function of the ground state of He-II is 
sought in the form^iiii 



(4) 



So = '^(^{k)p]^P--k+ ^ -= Pkl+k2P-klP-k2- 

k^^O ki.ksT^O ^ 

(5) 

The corrections of higher orders to Sq [Eq. ([5])] are ne- 
glected. In this approximation, the relation 



/(ki,k2) = - 



2cr(fci)2g(fc2)kik2 

e(ki-|-k2) + e(A:i)-|-e(fc2)' 



where 



e(k) = fc2(i _4fj(A:)), 



(6) 



(7) 



is valid^^. In works^"', an equation that connects 
[Eqs. (g]) and ([5])] with the He-II structure factor S{k) 
was derived. We shall write down this equation in ap- 
proximation ([5]) for 'I'o and using the notations of work^^ 
as follows: 



Aaiq) = 1 - 



1 



S{q) 



(8) 



where 



J_ ^ 8a(fc)a(k + q) + i?(k, q) 
^^„[l-Mfc)][l-4a(k + q)]' 



R{k, q) = 4/,(k, q) [1 + 2/,(k, q)] , (10) 



/,(k,q) -/(k,q) + /(-k-q,k)+/(-k-q,q). (11) 

Equations ([5]) and © were derived in work^"' from the 
known equation, which connects S{k) with the pair dis- 
tribution function F2{r): 



S{k) = l + n I (F2(r) - 1) e^^'^dr. 



(12) 



where n is the concentration of helium atoms. 

We note that in Eqs. ([8]) and ([9]), the interaction be- 
tween He^ atoms does not present explicitly, and 
is connected with the He-II structure factor only, so 
that such a way of finding ^'o allows the problem of 
descriptions of interaction between He* atoms at small 
distances^iiii to be avoided partially: provided strong 
overlapping of He* atoms, the description of atomic in- 
teraction using the interaction potential becomes inac- 
curate, because the atomic structure becomes impor- 
tant under such conditions, and it is necessary, generally 
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FIG. 1: Dependences of the quantity ka on k, k being in terms 
of fco — 2Tv/d — 1.756 units and d the average distance 
between He-II atoms. Squares mark a zeroth-order approxi- 
mation [Eq. p4|) ] for a-(k), pluses a first approximation, and 
triangles a second approximation calculated from the model 
"elliptic" potential with 17(0) = 60 Kii. 



speaking, to solve a quantum-mechanical problem of in- 
teraction of two nuclei and four electrons. The function 
*o [Eqs. (HI) and ©], found from Eqs. ®-([ni), takes 
short-range correlations into account more correctly than 
that found from the model potentialii. For a quite cor- 
rect account of the atomic structure, one should deter- 
mine the function ^'o for a system of nuclei and electrons 
rather than N structureless particles. It is a hopeless 
task. Nevertheless, as one can see below, configurations 
with atom overlapping are very improbable, so that from 
the physical point of view, it is quite reasonable to con- 
sider atoms as structureless particles. 

A single shortcoming made in the course of derivation 
of '^0 from Eqs. (|i))- pd|) was the break of series dH). But, 
since the model does not contain fitting parameters, the 
accuracy of approximation ([5]) can be estimated by com- 
paring both the theoretical spectrum of He-II quasiparti- 
cles and the theoretical potential of interaction between 
He* atoms with experimental ones. 

In order to find the wave function of the ground state, 
one has to know S{k) at the temperature T = K. As 
far as we know, the most exact measurements of S{k) 
were carried out in work2&. We used the smoothed data 
on S{k), obtained at T = 1 K in—-, and calculated the 
dependences S{k, T = 0) by the formula^^. 



S{k,T ^ 0) ^ S{k,T) tanh 



m 

2kBT- 



(13) 



At A: < 0.2 A-\ we supposed that S{k,T ^ 0) ^ k 
(because S{k = Q, T = 0) = O^^ and E{k ^ 0) = cfc in 
Eq. 1131)). In works^'^^, the integral equations ^ and 
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Q were not solved and a zeroth-order approximation 

1 



Aa{q) = 1 - 



S{q) 



(14) 



was used to determine ^'o- 
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FIG. 2: Function Si (r/a) [Eq. where a = 2.64 A is the 

"diameter"~ of the He* atom. The notations correspond to 
the same approximations for cr(fe) as in Fig. 1. 

Below, Eqs. ([6l)~(|ll|) will be solved numerically, and 
the solution cr(fc), which includes a single correction to 
the zeroth-order approximation (|14p . will be obtained; 
therefore, we shall call this solution a first approximation 
to CT(fc). The solution of the integral equation ([5]) cannot 
be found by the iteration method, so that we used the 
Newton one?^ for this purpose. As a result, two solu- 
tions were obtained, one of which, with a smaller energy 
per atom Eq — — 1.4 K (Eq = 0.1 K in a zeroth-order ap- 
proximation and —7.16 K in the experiment), being taken 
as the ground state. This solution for cr(fc) is shown in 
Fig. 1. 

A significant body of information concerning the prop- 
erties of \I'o is included into the function 



where 



k 



k#0 



(15) 



(16) 



The behavior of Si{r) at r — + shows how quickly the 
function decays if the atoms overlap. Fig. 2 repre- 
sents the function Si (r) for a zeroth-order (fT4|) and a first 
approximation for cr(fc), and for a{k) found in a second 
approximation, starting from the model potential^. One 
can see that 5i(0) ~ —1.7 in all those cases. One can sep- 
arate a two-particle summand of the form ^ o'(fc)pkP-k 

k#0 



from the addend in the r.h.s. of seeii, but a calcu- 
lation shows that the account of d-{k) renormalizes S'i(r) 
very slightly, by a few percent only. Thus, provided that 
two He'* atoms overlap, the wave function of the ground 
state diminishes by a factor of e^^'^ « 30, so that a sharp 
reduction does not occur, although the probability den- 
sity l^'oP decreases rather strongly, by a factor of 1000. It 
means that the He* atom possesses properties which are 
intermediate between "soft"- and "hard" -core ones. In 
case, for example, that 10 pairs of atoms overlap, ^'o de- 
creases by a factor of e^'^ ~ 10*'' as compared to its value 
for a uniform distribution of atoms without overlapping. 
Therefore, configurations where many atoms overlap are 
extremely improbable. 



III. CALCULATION OF THE HE-II 
QUASIPARTICLE SPECTRUM 

Knowing ^^o, one can find the wave function ^'k = 
V'k^o, which describes the state of the system with a 
single quasiparticle of the phonon type, and the quasi- 
particle spectrum E{k) from the following equationsi^: 



V'k 



ki#0,k , , 

P(k,ki) 

-p-k+ 2^ — pki-kp-ki+ 

ki ^ 

'--^^'- g(k,k„k,) 

TV 

ki,k2#0 



Pki+k2-k/5-kiP-k2 



(17) 



E{k)=e{k)+ J dkiP(k,ki)2ki(k-ki) + 

+ J dki(-2fc2)[g(k,ki,-ki) + 2g(k,k,ki)] 



(18) 



P(k,ki) [e(fci) + e(k-ki)-P(fc)] + y"dk2P(k,ki,k2) 
= 4a(fci)kki+2fc2/.(ki,k-ki), (19) 

F(k, ki, k2) = 4k2(k - ki - k2)Q(k, ki, k2)-f 

+ 2k2(ki-k2)g(k,ki-k2,k2), (20) 

g(k,ki,k2) [e(fci) + e(fc2)+e(k-ki-k2)-P(fc) = 
= P«(k, ki + k2) * G(ki, k2) + L(k, ki, k2), (21) 

F,(ki,k2) =P(ki,k2)-|-P(ki,ki-k2), (22) 

G(ki, k2) = [2a(fci)ki + 2a(fc2)k2] (ki + k2) + 

-^2(ki+k2)V.(ki,k2), (23) 
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FIG. 3: Theoretical He-II quasiparticle spectrum. Pluses cor- 
respond to a zeroth-order approximation (|25p for $o and ?/>k, 
triangles to a zeroth-order approximation for and a first 
one for i/^k, circles to a first approximation for and a second 
one for T/jk; the solid curve is drawn using the spline method; 
the dotted curve corresponds to the experimental spectrum^S,. 



L(k,ki,k2) -2k(ki+k2)/,(ki,k2). 



(24) 



In (|17p . the corrections of higher orders to T/jk are ne- 
glected. In Eqs. (fT8l) - ((24l we converted to the dimen- 
sionless variables fc' = fc/fco and Eik') — ^j^^r^, where 

fco = 2Tr/d and d ~ 3.578 A is the average interatomic 
distance. The primes will be omitted below. 

The quasiparticle spectrum, calculated using ^'o and 
-^k in a zeroth-order approximation (i.e. /(ki,k2) = 
0, Eq. (fT4|) . and i/'k = P-k), has the form of Feyn- 
man's known formulai, which describes the spectrum of 
a slightly nonideal Bose gas^^, 



21,2 



2mS{k)' 



(25) 



This spectrum is represented by pluses in Fig. 3. 

To find Tpk in a first approximation, we should assume 
that (5(k, ki,k2) = 0. From Eq. we have 



F(k,ki) = 



4a(fci)kki +2fcV,(ki,k~ki) 
e(fci)-f e(k-ki)-£:(fc) 



(26) 



The system of equations (IT51) and (pS)) was solved by the 
iteration method. The obtained quasiparticle spectrum, 
for cr(fc) in a zeroth-order approximation, is shown in 
Fig. 3. The spectrum of He-II in the indicated approxi- 
mations has been found earlier in work*. 

We note that the relation P(k, ki) ~ at k ^ and 
small ki is valid for ^p]^ in a first approximation, which 
corresponds to the Feynman-Cohen formula 



In a second approximation, it is necessary to solve 
the complete system of equations ([T ^ - (PH) . Similarly 
to work-ii, we solved these equations numerically. The 
system of equations ([T5|) . as a whole was solved by 
the iteration method, while Eq. (fTg]) by the method of 
quadratures^. In so doing, we used the values of cr^k) 
obtained in a first approximation. The error of the nu- 
merical definition of E{k) was about ±10 %; another 
error of about ±6 % stemmed from measuring S{k) in2& 
with an accuracy of ±2 %. The obtained spectrum E{k) 
is shown in Fig. 3. 

/^From Fig. 3, one can see that if the number of correc- 
tions, which are taken into account, increases, the agree- 
ment between the theoretical and experimental spectra 
improves, so that for and ipk determined in a first 
and a second approximation, respectively, we have a 
good agreement between E{k) and the experiment. The 
"shoulder" E{k) « 17 K in the experimental spectrum at 
k > 2.5 A is connected, in our opinion, with a hybridiza- 
tion of the spectrum that describes a single quasiparticle 
with a two-roton level-^. 

Knowing the structure factor, one can restore the inter- 
action potential between He'* atoms by finding cr(fc) and 
/(ki, k2) from Eqs. ([S|)" ([TT|) with known S{k) and substi- 
tuting the obtained solutions into the following equation 
for the Fourier-image ^{k) of the potentiaUi: 



-<7(ki)kl -t- 



ni'{ki 



<y\k,)kl 



_ 1 
~ N 

The potential 



k2#0,-ki 



/,(ki,k2)0.5(A:| 



kik2 



U(r) 



(27r)3 



iy{k)e'^''dk, 



(27) 



(28) 



where i/(fc) is a solution of Eq. ([?f|) . is shown in Fig. 4 
for o'(fc) taken in a zeroth-order and a first approxima- 
tion. The potential U (r) in a zeroth-order approxima- 
tion was obtained in work^^ earlier. The potential cal- 
culated by us approximately agrees with those obtained 
in worka^iiii^i^i, but not with Aziz's potential^, which 
possesses a very high barrier of repulsion J7(r = 0) ~ 
10^ K. This discrepancy might be caused by the ef- 
ficiency of the potential that describes the interaction 
between He"* atoms at small distances, as well as by dif- 
ferent modeling of such interaction. It is not improbable 
that some processes (e.g., the scattering of He* atoms) 
are better described by Aziz's potential, while others (in 
particular, the calculation of '^q, and the E(k) spec- 
trum) by a potential with a much smaller effective barrier 
U{0) ^ 100 K. One can see from Fig. 4 that the found 
potential has a minimum at rmin = 3 A with the depth 
Umin = — 7.7 K, which approximately corresponds to the 
Lennard-Jones experimental "well" with rmin = 2.97 A 
and [/,„i„ = -10.8 K35. 

Fig. 4 also testifies that U{r) for a{k) in a first ap- 
proximation differs appreciably from that for a'(fc) in a 
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FIG. 4: Potential of interaction U{r) between He* atoms, re- 
stored from the structure factor. The sohd curve corresponds 
to a zeroth-order approximation for (T(fc), and the circles to a 
first one. 



FIG. 5: Fourier image V(k) = ni'{k) for potentials U{r), 
shown in Fig. 4, ko = 2n/d. The notations correspond to the 
same approximations for a{k) as in Fig. 4. 



zeroth-order one. The potential changes appreciably even 
if Eqs. ^ and (O are rewritten in another but equivalent 
form. The inaccuracy of the U (r) determination stems 
from the fact that, according to Eqs. (^7)) and (|28p . the 
potential U{r) depends strongly on the a{k) values at fc's 
that are not small, k = 2fco 4fco, because the summand 
J a{k)k'^dk makes a contribution to U{r). The values of 
cr(fc) are small at such fc's, but the corrections to (j{k) 
turn out to be of about cr(fc) itself. Therefore, (T(fc) is 
not determined exactly at considerable fc; this circum- 
stance has almost no influence on the resulting quasi- 
particle spectrum, but induces a significant error while 
finding U{r). Thus, we can only estimate the potential 
U{r), but in order to calculate U{r) with a higher ac- 
curacy, one must determine the next approximations for 
cr(fc) and measure 5'(fc) more precisely. 

IV. COMPARISON OF DIFFERENT HE-II 
MODELS 

Below, we present a short, schematic comparison of 
various methods which are applied in order to explain the 
microstructure of He-II. In so doing, we do not pretend 
that our analysis is complete or perfect. 

There are plenty of works dealing with the micro- 
scopic description of He-II. Some analysis can be found 
in reviewsi^ii^. The main approaches are as follows. 

(i) Semi-phcnomenological methods, where certain 
equations (like the Gross-Pitayevskii one^l or that 
of the model of a "continuous medium"—) are pos- 
tulated and used as a start point to derive the 
quasiparticle spectrum. Several fitting parameters 
(FPs) are used at that. The main shortcoming of 



these methods is that it is not clear how precisely 
the initial postulates correspond to the He-II mi- 
crostructure. 

(ii) Microscopical approaches, which are based on the 
calculation of ^'o and V'k- 

(a) A "straightforward" solution of the A^-particle 
Schrodinger equation in the or the k-space 
(the MCV, see Refs.^'^'iajJ. and this work). 
Here, both and V'k can be determined with- 
out introducing FPs. 

(b) "Indirect" methods for solving the 
Schrodinger equation, e.g., the "corre- 
lated basis function"— and "hypernetted 
chain"— 1^ approaches. 

(c) Variational methodsi'^'^i^'^. 

(d) In the "shadow wave function" (SWF) 
approacbi2ii2iii^iS£, the attempt is made to 
partially "contract" the whole infinite series 
of correlative corrections to In and ^'k into 
separate simple "shadow" factors. This pro- 
cedure was argued by a certain reasoning, in 
particular, by that taking the delocalization of 
atoms into account means the partial consid- 
eration of higher correlations. As was noted 
inii, such a solution is the first iteration of 
the Schrodinger equation represented in the 
form of a functional integral. Drawbacks of 
the models are as follows: the exact solutions 
for ln\l/o and 4'k are infinite series, and it is 
not clear to a which extent the shadow factors 
will allow one to evaluate this series; moreover, 
too much FPs are in use. 
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(e) The numerical Monte-Carlo (MC) 
method^'— which gives the most ex- 
act description of the ground state, its energy 
Eq, the structure factor S{k), and the values 
of all condensates. But the method does not 
allow one to see the analytic structure of a 
solution and its details and does not yield the 
curve E{k). 

In approaches c, d, several FPs are used. The main 
lack of all models a-d from (II) consists in that the 
exact solutions for In and ■(/'k are infinite corre- 
lation series. In practice, one succeeds to consider 
only 2-3 first terms, whereas the omitted correc- 
tions are not small. 

(iii) Field-theoretic models. 

(a) Studies of the total Hamiltonian H in the k- 
space^ii^iiiii^i^i^. To a certain extent, this 
case is rather close to item (ii,a). The conden- 
sates do not appear explicitly in the equations. 

(b) Studies of the Hamiltonian H in the k-space, 
in the representation of the operators and 
Sk for quasiparticles^^i^i^. Here, the conden- 
sates appear explicitly. 

(c) Solution of the equations similar to the 
Belyaev-Dyson ones^i^. 

Models b and c involve FPs. 



In our opinion, the most perspective may be the field- 
theoretic approaches of types (iii,b) and (iii,c), the MC 
method, or quantum-mechanical methods that have not 
been discovered yet, which will start from exact micro- 
scopic equations, will not use fitting parameters, and 
where the expansion in a small parameter will be carried 
out. At the same time, approaches (ii) form a necessary 
complement to (iii) ones. 



V. CONCLUSIONS 

To summarize, in this work, using the method of col- 
lective variables, the spectrum of He-II quasiparticles has 
been obtained, and the wave functions of the ground and 
a first excited state of helium-II have been found approx- 
imately, without introducing any fitting parameter into 
the model. We have solved the equations that had been 
derived from the exact microscopic equations. A single 
inaccuracy of the method consisted in breaking the series 
for and ^Ak- The obtained He-II quasiparticle spec- 
trum agrees well with the experimental one. Therefore, 
we believe that the found solution reflects the microstruc- 
ture of He-II. This result makes the solution found in 
work^ earlier more accurate. 

The author is grateful to V.E. Kireev for discussion of 
numerical methods and to E.A. Pashitskii for discussion 
of the work and useful criticism of previous results. 
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